Ab-initio formulation of the 4-point conductance of interacting electronic systems 
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We derive an expression for the 4-point conductance of a general quantum junction in terms of the density 
response function. Our formulation allows us to show that the 4-point conductance of an interacting electronic 
system possessing either a geometrical constriction and/or an opaque barrier becomes identical to the macro- 
scopically measurable 2-point conductance. Within time-dependent density-functional theory the formulation 
leads to a direct identification of the functional form of the exchange-correlation kernel that is important for the 
conductance. We demonstrate the practical implementation of our formula for a metal-vacuum-metal interface. 
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I. INTRODUCTION 

Impressive progress has been achieved within the non- 
equilibrium Green's function (NEGF) formulation of quan- 
tum transport using the simple ground-state density-functional 
exchange-correlation potential in a self-consistent formula- 
tion'-^ (NEGF-DFT). However, limitations of the latter ap- 
proximation were recently identifiedMi^j^iLL2ii2. For in- 
stance, NEGF-DFT's omission of the derivative-discontinuity 
in the exchange-correlation energy functional was found re- 
sponsible for serious errors in transport calculation through 
localized resonant levels^ *'. Improvements through an (spin- 
) unrestricted NEGF-DFT formulation have been argued to 
describe properly some aspects of the Coulomb blockade in 
quantum junctions^. At this level of the theory the exchange- 
correlation potential of the equilibrium system, v^^c, is respon- 
sible for the electron interaction effects. 

In a further theoretical development, Na Sai et al.^ iden- 
tified a dynamical correction to the resistance of a quantum 
junction stemming from the contribution of the exchange- 
correlation electric field to the overall drop in the total po- 
tential. They estimated the coiTection within time-dependent 
cuiTent-density functional theory" (TDCDFT) and showed 
that it has its origin in the non-local density-dependence of the 
functional. The very applicability of time-dependent density- 
functional theory (TDDFT) to the problem of quantum trans- 
port in the long-time limit has been discussed in depth by Ste- 
fanucci and Almbladh^ and by Di Ventra and Todorov^. 

Several authors have proposed alternative treatments that 
avoid the complexities of the exchange-coiTelation kernels 
of TD(C)DFT, either by treating the central region with the 
configuration integration me thod'^''^ while approximating the 
non-equilibrium distribution of the electrons, or by using the 
usual NEGF-DFT approach in combination with a model self- 
energy within the central region'**. More systematic approach 
to the self-energy can be obtained through well tested approx- 
imations like the GW method '^■'^. Due to the large computa- 
tional demand so far only very small systems with restricted 
size of the basis set could be studied. However, the results are 
encouraging, e.g. the Kondo effect phenomenology seems to 



be well described within the self-consistent GW methodic. 

Nonetheless, a systematic approach for addressing the con- 
ductance of a fully interacting system at the ab-initio level 
is not available. This is partly due to the fact that the very 
formulation of the NEGF formalism'^ is based on the con- 
cept of noninteracting electrodes and demands partitioning of 
the system and the electron-electron interaction'**. Stefanucci 
and Almbladh^ showed that the partitioning can be avoided in 
principle but practical inclusion of the many -body interactions 
into the formalism seems to be extremely cumbersome. Fortu- 
nately, partitioning is not necessary within the linear response 
formulation. Several authors have addressed the conductance 
of interacting system of electrons within the framework of the 
Kubo formalism arriving at a 2-point Landauer-like formula 
for the conductance by making certain assumptions 

about the steady-state total electric field. There is, however, 
a problem with this derivation since it ignores the charge re- 
distribution in the conducting system when the steady state is 
forming. In fact, ignoring these aspects one quickly arrives 
at various unphysical coiTections to the conductance"^. The 
problem of charge redistributions has been first pointed out 
by Thouless^. Later Kamenev and Kohn^** cast it into a self- 
consistent framework for many- terminal conductance. 

In our work we further develop the formalism that treats the 
charge redistribution correctly, find its physical interpretation 
in terms of a time-dependent transient process that leads to 
the establishment of a cuiTent-caiTying steady state, and de- 
rive a closed formula for the 4-point conductance that is well- 
defined formally as well as physically. 

To describe quantum transport it is necessary to consider 
the production of a steady state with nonzero current in a 
generic quantum junction, starting from the ground state (or 
more generally, from an equilibrium ensemble at finite tem- 
perature). Physically, we expect to achieve a steady state that 
has a uniform average current / flowing through the system, 
accompanied by a total potential exhibiting a drop V across 
it. As was advocated many years ago by Landaue?^ V arises 
from the resistivity dipole, a local charge imbalance of a piled- 
up charge in the front and depleted charge behind a quantum 
junction. Physically such a state is reached by attaching a 
macroscopic battery to the circuit that pumps electrons from 



one electrode to the other until a certain potential difference 
between these two is established. For modeling purposes this 
mechanism cannot be used directly^. Here we will use a 
construct based on an auxiliary homogeneous vector poten- 
tial A""''(f ) = - /o'£"'"(f')c/f', pointing in the direction of the 
eventual flow of current, giving a momentum transfer to all 
the electrons in the infinite system within a finite interval of 
time. (This auxiliary vector potential is similar in spirit to the 
one used by R. Gebauer and R. Car to model a nanojunction 
using periodic boundary conditions'^. There, in contrast to 
our work, the vector potential grows linearly in time for all 
f > and the work exerted on the system must be dissipated 
via auxiliary phonons located in the electrodes.) The infinite 
extent of the system is formally essential to our treatment. It 
is necessary for a continuous spectrum and for giving a mo- 
mentum transfer at f = to infinitely many electrons present 
in the system so that the current will flow for all times t > 0. 
However, we would like to mention already at this point that 
in practical calculations the necessity of the infinite extension 
can be relaxed so that the formulation can be used for practical 
ab-initio calculations. 




FIG. 1: (Color-online) Measurement of the voltage between the 
macroscopic electrodes, where the current density is zero, gives the 
two-point conductance G^^ = // (A/i /e) . Measuring the voltage drop 
inside the simulation box (supercell), where the current density is 
nonzero, gives the four-point conductance G"*^ = I/V. Increasing 
the supercell the two quantities approach each other (see SeclVTl. 



II. THE NON-LOCAL CONDUCTIVITY AND 
CONDUCTANCE 

The response of the current density to a general external 
electric field is given within the linear response theory by'^ 

j{r,t)^J' dt' J d^r'd(f,r';t-t')-E''^(f'y). (1) 

For simplicity we introduce a symbolic notation for the above 
equation in the form / = a*E^^. In our work, the role of the 
external field is taken by the auxiliary field that is homo- 
geneous and has only a finite duration, i.e. it is absent for large 
times and it bears no information about the drop in potential in 
the long-time limit. The latter is contained within the induced 
field E'{r,t), which appears explicitly when considering the 
irreducible (or proper) conductivit y'^''^ , 

j = d* E""' = a'"' * + £'■)• (2) 

The induced field accounts for the electron-electron interac- 
tion at the Hartree level, which is usually referred to as the 
long-range effects of the Coulomb interaction, whereas the 
rest of the interactions between electrons, the short range part, 
is included within the irreducible conductivity. The quantity 
in which we are primarily interested in is the conductance, 
defined as G = I/V where V is a voltage drop. For G, the de- 
tailed spatial and time dependence of the current density and 
the induced field while the steady current is being established 
are of no relevance. We define the voltage drop as the overall 
drop in potential of the induced field E' along the current flow 
(along the z axis) between far left and far right, 

y = lim dzE[{r,t), (3) 



which will be indicated by a superscript 4P in the associated 
4-point conductance, G^^ . This definition is most suitable 
for ab-initio modeling and does not suffer from ambiguities 
present when one considers two chemical potentials. In meso- 
scopic physics, the above definition of conductance is referred 
to as the 4-point conductance (hence the superscript) since 
it represents measurement of the voltage drop using contacts 
different from those acting as a source and drain of the cur- 
rent. In the extreme case of a ID conducting channel of non- 
interacting, but locally neutral electrons it reduces to the well- 
known expression — {2e^ / h){T / R). This is in contrast 
with the most frequently encountered 2-point conductance, 
G'^, where the voltage drop is understood to be the difference 
in electro-chemical potentials of the two macroscopic elec- 
trodes, A/i/e (see Fig.[T]). Since the electrodes are not part of 
the quantum-mechanical model of the conducting channel, the 
familiar Landauer result G'^ = {2e^ /h)T is believed not to be 
derivable from the Kubo formalism^" and its plausible deriva- 
tions are always accompanied with steps motivated by phys- 
ical insight and arguments about phase-incoherent adiabati- 
cally widening electrodes'^^. Here we show that for a nanocon- 
tact between massive electrodes the 4-point conductance ap- 
proaches the 2-point Landauer formula. However we never 
make a reference to A^u/e and we consistently work with the 
drop in the induced electrostatic potential only. Our claims for 
its applicability for the experimental 2-point conductance fol- 
low from the fact that for massive electrodes connected with a 
nanojunction the current density in the electrodes goes to zero 
and thereby the drop in the induced (Hartree) potential is an 
excellent indicator of the phenomenological quantity A/i/e. 
The latter argument serves as a motivation to use the G'*^ as 
conceptually and formally well defined quantity to character- 
ize systems of fully interacting electrons in any geometry of 
the junction. 
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III. SINGULARITY IN THE RESPONSE FUNCTION AND 
THE 2-POINT CONDUCTANCE 

An essential property of the conductivity response function 
of extended systems capable of coherent transport is their sin- 
gular long range character for large times or - equivalently - 
as the frequency approaches zerc^' . To show this behavior ex- 
plicitly we need to introduce a few geometrical details of the 
studied nanojunction which do not restrict the generality of 
our argument. We consider a geometry where z is a direction 
of the current flow, q a reciprocal vector in that direction, and 
for clarity we assume that the system can be put into a peri- 
odic supercell along the x^y directions. The conductance G^^ 
is then a conductance for the area of the supercell S. To find 
the current flowing through S, we integrate Eq. Q across the 
cross-sectional area. This naturally leads to the cross-section 
integrated conductivity ^*^ 

o'"{z,^-t)^ I [ dS-&"'{r,r'-t)-dS', (4) 
JsJs 

which relates the current to the z-component of the total field, 

I{z,t) = fdt'J°° dz'a"-'-{z,z'-J-t'){E"'"it')+E\z',t')). 

(5) 

The singular long-range character appears as the indepen- 
dence of G"'''{z,z'',t) on z and z' as f ^ oo which when per- 
forming the Fourier transforms z,z' q,q' and t ^ co + ia, 
takes the form^^ 

lim a'''''{q,q',co + ia) = 27i:G^^5{q)5{q'), (6) 

where we have formally introduced the 2-point conductance, 
G^^, as the strength of the above mentioned singular character. 
To extract the strength we introduce a linear functional [] 

^"[a"-^] = lim / [^a'"-{q,q'Ua), (7) 

for which we simply have G^^ = ^^'^[ct'"']. It has been shown 
previously that this identification is correct for noninteracting 
electrons^'. Further motivation to refer to it as a 2-point con- 
ductance for interacting systems will be discussed in Section 

lYl 

IV. ONSET OF THE CURRENT INSIDE THE 
ELECTRODES 

The Eq. (|5]l can now be analyzed in the light of the above 
observations. First, let us consider a distant part of one of the 
electrodes, |z| 3> 0, for times such that XjEp -C f <C |z|/vf , 
where Ep and vp are the Fermi energy and Fermi speed re- 
spectively. In this regime, the response in the electrodes 
(|z| ^ 0) is determined by the electrode itself, independently 
of the nanojunction. After a short relaxation time (~ 
the auxiliary field establishes a uniform (i.e. z-independent) 
current inside the electrode 

/(z,f) = /'^z'/' dt'&"\z,£;t-t')E""\t'), (8) 



where &^^''' is the irreducible conductivity of the electrode 
only'*^ while preserving the local charge neutrality within the 
electrode. We are free to choose any time-dependence of the 
field as long as it delivers finite change of momentum for each 
electron and thus establishes the eventual steady current; the 
most convenient choice is E"'"{t) = —A""^8 (f ), where A""'' is 
the magnitude of the change in the homogeneous vector po- 
tential at f = 0. 

This picture is not valid as one approaches the region of 
the nanojunction, i.e. for z ^ 0. Still, since the current will 
keep coming from the left electrode and disappearing into the 
right electrode (there is always such |z| ^ that for any time 
t we have t ^ \z\/vf and hence the above analysis can be 
applied), the charge will have to pile up in front of the junction 
and similarly deplete behind the junction, creating a resistivity 
dipole^^. The incoming current cannot be decreased by some 
form of reflected front/disturbance from the junction since this 
would break the local charge neutrality within the electrode, 
leading to strong opposing fields. Hence the dipole around 
the junction and therefore the drop in the induced potential V 
will grow, increasing the current in the junction region until 
it becomes equal to the current, /, deep inside the electrode, 
as described in the preceding paragraph. This will be possible 
since in the linear regime we expect / — G'^^V. 

Having in mind this physical process, the established uni- 
form current at long times expressed by Eq. ([8]l for the distant 
electrode applies at all positions including z — 0. We express 
the electrode conductivity cj"^'^''' in reciprocal space for long 
times as 

/(f^oc) = - lim f dqa"'''''(q,q' = 0;ia)A''"' (9) 

= ^'^[(7''"''-"]A"", (10) 

where we use a small imaginary frequency ia to perform the 
long-time limit^i. The linear functional [] will be further 
discussed in Section lVIll 



V. THE 4-POINT FORMULATION OF THE 
CONDUCTANCE 

The above physical picture can be directly used within 
Eq. Q. We write the conductivity of the whole system as 
the conductivity of an infinitely long electrode alone, a"'''''', 
plus a further term, cr'"^^, characterizing the presence of the 
junction: g"^'^ = o'^''^'^ + ff'" -'. Eq. (|5]l is thereby cast into the 
form 

/ ^ (7''''"'^*£'"" + (7''"'''^*£"'"+(7'''"*£'. (11) 

Evaluating this equation at z — and using the definitions of 
the functionals ■^'^ and the equation (fTTT i can be written 

as 

/ = ^'^[a''-'"'«]A™-'^ + ^'^[a''-'"J]A"'">^^[a'"']V. (12) 

From Eq. (|8]l we know that in the long-time limit, the first 
term on the RHS is itself equal to /, so that the second and 
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third terms sum to zero, from which we easily arrive at the 
central result of this paper 



g4^ = 



(13) 



where we have used the fact that [a"''-''] = ^°"[cr"''] - 
^cT[(y,>r,e] jjjjj jjj^j y ^ JdzE'iz) = = 0). 

It is very instructive to apply the above result for a simple 
ID non-interacting gas with a single scattering center'*'*. For 
this system, the noninteracting conductivity in the small fre- 
quency limit has a form^'* 



a^{q,q';ia) = 



Vftt 



2 

-R- 



f a 



:5{q-q') 



Kv^q^ + a^ nv],{q'Y- 



(14) 



The irreducible conductivity of the homogeneous ID gas, cor- 
responding here to the electrode in the general case is simply 
the first part of the above expression 



a ''^{q,q';M) = — - 



a 



:5{q-q'). 



(15) 



It is straightforward to see that using the above forms we 
have ^'"[ff'^] = (1 -R)/n = T/n, .^""[a^''] - vp/jta and 
^cjfjOjj ^ Rvf/na, and we arrive at the 4-point conduc- 
tance G'*^ — {1/%){T /R). The arguments of local charge neu- 
trality in the electrodes and its consequences are very closely 
related to the original comment by Thouless^- and a more re- 
cent work by Kamenev and Kohn^'*. In the form presented in 
Eq. ( fT3] l it is valid for nanojunctions of any shape as long as 
we can put them into a supercell, including the extreme case 
of a planar metal-vacuum-metal junction for which we will 
demonstrate its applicability in Section lVlIl 

The generality of Eq. ( fT3] l can be used to show the equiva- 
lence of the 2-point and 4-point conductances for nanojunc- 
tions with massive electrodes. The demonstration is based 
on the observation that for such a geometry [&"'■"] ^ 
,j^"[o'"'] and hence the prefactor containing the =^[]s in 
Eq. ( fT3] l goes to 1. The mentioned inequality is immedi- 
ately clear from the following argument: [ff '"'■''] gives the 
current flowing in the electrodes when the vector potential is 
changed from zero to A"'", whereas [f'" ] formally gives 
current as a response to the same disturbance but in a fictitious 
system having the same geometry as the real nanoj unction but 
for which the long-range Coulomb interaction is missing. This 
implies that the local charge neutrality is not enforced (since 
only the irreducible response enters ^°^) and as a result most 
of the electrons coming from one electrode will be reflected 
from the nanojunction and decrease the current. Hence the re- 
sulting current will be much smaller in this second case and 
the inequality is fulfilled. This formally justifies Landauer's 
arguments in favor of the expression G~T/nas the conduc- 
tance of a nanojunction between two phase-randomizing and 
adiabatically widening reservoirs. As we can clearly see, it is 
only the widening that is really needed to have the conduc- 
tance of the junction as a whole be given by Eq. Q. This ex- 
plains our motivation to refer to the latter quantity as a 2-point 



conductance even for interacting electronic systems. The ar- 
gument applies also to opaque barriers, not necessarily having 
a constriction, which has been explored in previous works^^*^ 
and is also demonstrated in Section lVlIl 



VI. INCLUSION OF EXCHANGE AND CORRELATION 
INTO CONDUCTANCE 

Having established the validity and generality of Eq. (fTSl l 
we can proceed to approximations that go beyond non- 
interacting or Hartree-like interacting electrons. The sim- 
plest way to do so is within the framework of time-dependent 
density-functional theory. Some doubt about the applicability 
of the latter might arise from the dynamics deep in the elec- 
trode described in Section |IV] There, the induced current is 
dominated by a divergenceless component that cannot be re- 
lated to a time-dependent density, which is the only physically 
relevant quantity within the TDDFT. However, at the level of 
Fermi-liquid theory, this divergenceless current is identical to 
that of a non-interacting system due to the backflow of quasi- 
particles^^. Hence for this part we do not expect corrections 
arising from the exchange and correlation and we should ex- 
pect ^^[a"''''^] = ^'^[a^^-^], where a^'" is the non-interacting 
conductivity of the electrode. On the other hand, the response 
described by the rest of Eq. (fTTT i. leading to the voltage drop in 
the induced potential, is essentially localized around the nano- 
junction and is completely described by the time- and space- 
dependence of the electronic density which is amenable to the 
TDDFT approach. 

To cast our theory into the TDDFT framework, we need to 
reformulate the functionals and '^'^ since the conductivity 
is not directly accessible within TDDFT. The irreducible con- 
ductivity is simply related to the irreducible polarization^ij^^ 



a 



a"'''{q,q',ia) = -x"'''{q,q';ia), 

qq' 



(16) 



which is conveniently calculated via the density response 
function calculation within TDDFT. To incorporate this we 
also introduce new functionals of the irreducible polarization 



a 



qq 



a 



^^"[ X'{^,q'-Ja)]. 

qq' 



(17) 
(18) 



for which we give explicit expressions suitable for direct nu- 
merical evaluation in reciprocal and real space representation 
in the Appendix. The irreducible polarization satisfies the 
Dyson equation 



X"'''{q,q'\ia)=X^{q,q'\ia) + 
dq"dq'"x' (q, q"; ia)Uc {q" , q" \ ia)x'"' {q'\q'; /a) , (19) 



where X^{l^l''->^(^) is the non-interacting Kohn-Sham den- 
sity response function, and f^dq" ,q"'',icc) is a non-local 
frequency-dependent exchange-correlation kernel^''. The 
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above two equations give the 4-point conductance of an in- 
teracting electronic system: for a given kernel one needs 
to invert the Dyson equation (l9[ . substitute the result into 
Eq. ( fT6b and employ the general expression Eq. ( fTsT l. 

However, we can gain more insight by multiplying Eq. (fT9] l 
by —a/{qq') and taking the limit a ^ 0. Clearly, the re- 
sulting left-hand-side is singular in q, with the strength be- 
ing G^^ according to Eq. (|6]l. The strength of the first term 
on the right-hand-side of ( fT9t multiplied by the same factor 
gives the 2-point conductance of the non-interacting Kohn- 
Sham system, G^^". The difference between these the two, 
i.e. the correction due to the exchange-correlation kernel, 
is then nonzero only if the last term in Eq. M% also leads 
to a singular form. This observation can be used to de- 
duce the forms of the kernel that do influence the conduc- 
tance, since the general character of % for small a is well 
known. The most interesting choice, making use of the char- 
acter of x"/''' - qq', is f:,ciq,q';icc) = ■^A{q,q'\ia) where A 
is well-behaved: A{q,q';ia) R'^^" = A{q = 0,q' = 0)^0 
for a — > 0. The resulting 4-point conductance then takes the 
form 



(1 +'^[x°]R'^^).^[x^'''] - ■'P\X^\ 



\X\ (20) 



In the case of a narrow junction between massive electrodes 
and using the arguments leading to Eq. ( fTSl ), we obtain 



(21) 



i.e. in this very important case 7?'^^° represents part of the 
dynamical resistance, having origin in exchange-correlation 
effects and can be accounted for by adding this resistance in 
series with the Kohn-Sham result. 

An approximate correction of this form has, in fact, been 
identified by Na Sai et al? for interacting electronic systems 
with weakly inhomogeneous potential along the direction of 
the current flow. For such systems it is possible to show that a 
purely longitudinal exchange-correlation electric field used in 
their treatment within TDCDFT-ii^ is equivalent to a contri- 
bution to the TDDFT kernel of the asymptotic form for small 
q,ci 



m 



qq 



dz- 



4ri ( d.n{z) 



qq 



where A is the cross-sectional area of the considered system, 
n{z) is the electronic density and r\ is the dynamical viscos- 
ity of a homogeneous electron gas''. We should note that 
for homogeneous systems Z?"^^" = since d^n{z) = 0. This 

is important since the functional form f^c given above and 
the limiting process would not lead to a finite result for a ho- 
mogeneous system. It is interesting to note that the nonlocal 
character of the kernel, i.e. f^'^ ^ 1 /q^ has been found also re- 
sponsible for significant improvement in the ab-initio studies 
of the optical spectra in many materials^-^i^i. 



A very popular approximation for the exchange-correlation 
kernel is the adiabatic local-density approximation (ALDA)^ 



f^r'\ry)=f^r[n\vmv-v') 
fi^''^[n']=d\n'e.M'))ldn\ 



(23) 
(24) 



(Here evt («*') is the exchange-correlation energy per par- 
ticle of a homogeneous electron gas (HEG) of density 
«"), leading to a kernel in the reciprocal space of the 

form f^c^°^\q,q';(0) = fxc^°^\q - q'\(0) B{(o)5{q - 
q'),B{co) B ^ Q for o ^ 0. It has been argued by sev- 
eral authors that ALDA should not contribute to any change 
in the conductance.^^ The basis of this argument comes from 
the fact that one can account for the exchange and correla- 
tion effects either via the exchange correlation kernel, or by 
using the exchange-correlation electric field. The latter can 
be shown to contribute to the conductance only if there is a 
nonzero drop of the exchange-correlation potential across the 
system, which is absent within any local or semi-local approx- 
imations. The argument is valid also within our theory when 
applied to Eq. (fTTT i. However, it is interesting to explore the 
consequences if one stays within the kernel-based treatment. 
Here the proper definition of the 4-point conductance is es- 
sential. If one directly uses the 2-point formula in Eq. (|7]i for 
metal-vacuum-metal interfaces, one finds nontrivial ALDA 
corrections solely due to the fact that the Kohn-Sham conduc- 
tivity does redistribute the charge at long distances when a lo- 
calized Kohn-Sham potential with nonzero drop is applied^-. 
The numerical results for metal-vacuum-metal presented in 
Section lVlIl show a cancellation between these ALDA correc- 
tions within ] and those entering ,'^[x"''] within the pre- 
cision of the numerics so that the G^'^ coming from Eq. (fTsT l 
remains unaffected by the ALDA kernel, as it should be. It is 
important to emphasize this point because studies of exchange 
and correlation effects using the 2-point formula only for ID 
atomic chains^" or quantum wires^^ may easily lead to incor- 
rect conclusions. 



VII. CONDUCTANCE OF A METAL- VACUUM-METAL 
INTERFACE 

The metal-vacuum-metal junction is one of the simplest 
systems for ab-initio study of electronic transport and hence it 
conveniently serves as a demonstration that the formula ( fT3] l 
can be numerically implemented even for realistic systems. 
In practical implementations we cannot work with an infinite 
system, for which the response function that enters the func- 
tionals ] and ] is needed. Instead, we take a finite sys- 
tem with z G {—L,L) with zero boundary conditions at the end 
points. This brings certain restriction on the zero-frequency 
extrapolation that we will discuss below. The two perpendic- 
ular directions can be dealt with easily in the reciprocal space 
for this particular system. 

In our calculations we employ two jellium slabs of thick- 
ness I and density given by rj = 3.0 a.u., separated by a dis- 
tance d, i.e. 21 + d < 2L where the supercell extends from —L 
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FIG. 2: (Color-online) Dependence of G^''{a) ='Sa{T'. 
quency oc for several electrode widths / (vacuum is fixed for d = 
8 a.u.). For frequencies below the thick line, given by «„„■„ (/) = vp jl, 
the finite size effects appear, so that a ~ 0.05 are used for extrapola- 
tion to zero. The vertical shifts of the individual curves is caused by 
the oscillation of the Fermi energy (see also Fig. [3}. 
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FIG. 3: (Color-online) Dependence of the extrapolated G^^ on the 
electrodes' widths / for d = 8.0 a.u for a square-barrier potential. 
The inset shows the oscillations in the Fermi energy and the conduc- 
tance G^^ in a restricted range of / for a self-consistent calculation. 
Choosing widths for which the Fermi energy attains its minima (in- 
dicated with dots) leads to a stable estimate of the conductance of the 
infinite system. 



to L. The calculation of [z^tI \ia) is performed at the self- 
consistent LDA level'*'^. Subsequently for the ALDA calcula- 
tion we invert the Dyson equation, Eq. dUi, in real space and 
thereby calculate the irreducible response X^"' [z^t! \i(X) . The 
inclusion of the non-local kernel, such as that arising from 
the viscosity of the electron liquid, can be done directly via 
Eq. ( |20] |. avoiding the inversion of the Dyson equation. For 
the evaluation of ^[ ] and ] we employ their real space 
form (Appendix), integrating over the simulation cell, 

^^ate'H = « r ['^dzdz'x"'''{z,z';ia), (25) 
J-lJq 

^a[x'"] = afdz dz!x'"{z,z)z, (26) 

J-L J-L 

keeping a finite for the moment. For extrapolation we need 
to use imaginary frequencies a <C fif so that the transient dy- 
namics is removed from the response function. On the other 
hand, the finite extent of the electrodes restricts the limiting 
procedure for a ^ 0; one can go down no more than to imag- 
inary frequencies o:,„,„ « Vf // (vf « 0.64 a.u. for = 3 a.u.). 

To demonstrate the precision and scaling of our calculation 
we first present calculations for a non-selfconsistent square- 
barrier potential (width 8 a.u., height 0.25 a.u.) between 
two 3D electrodes (Ep = 0.2 a.u.). The chosen values are 
reasonably close to our self-consistent potential but at the 
same time allow for comparison of with the exact value of 
G^^ — lima^o [z"] using the analytically known form of 
the transmission probability. The calculated values of 'i^aix^] 
shown in the Fig. |2] for finite slabs and finite frequencies 
clearly show finite size effects below a„„„ . However, even for 
frequencies larger than a,„i„ we observe a remanent horizontal 
oscillatory dependence of the conductance curves on the elec- 
trodes' width. Even though the amplitude of the oscillations 
for smaller widths are substantial, the overall convergence to 
the exact value is evident. The oscillatory character of the 
Fermi energy can be effectively used for identifying the most 
suitable system widths for an optimal estimate of the infinite 



system conductance. One should choose the widths for which 
the Fermi energy attains local minima (inset in Fig.O, which 
is the method we use within our work. 

Similar extrapolation to zero frequency can be done for 
^[ ], or even better for ax ^[] which approaches a constant 
value for an infinite system. However, since our aim is to cal- 
culate the 4-point conductances, we directly extrapolate the 
expression (l20l i which behaves very similarly as the G^^(a) 
described in detail in the preceding paragraph. 

We now apply our methodology to a fully self-consistent 
calculation for a sequence of vacuum widths d ~ 1 , . . . , 10 a.u. 
to study the effect of the recently suggested viscosity-related 
exchange-correlation coiTection due to Na Sai et al?, and to 
explore the cancellation of ALDA corrections between ] 
and ]. The overall dependence of the conductance on d 
is exponential. For clarity of presentation we show first the 
exponentially decreasing form of the calculated Kohn-Sham 
conductance, Gq^, in the upper panel of Fig. HI the small sym- 
bols indicate error bars in Gq' due to uncertainty in the extrap- 
olation to zero frequency. For this system the absolute values 
of the corrected conductances after inclusion of either TDDFT 
kernel - ALDA kernel or viscosity kernel - lie within these 
indicated eiTor bars. As we have discussed in Section lVTl the 
ALDA does not lead to any correction, but numerically this re- 
sult is not trivial. In fact both, the ^^[z'*^™] and ^[j'*^™] do 
change due to the presence of the exchange-coiTelation kernel 
(l24l i. but these changes are canceled in the total expression for 
the 4-point conductance (fT3] l within the numerical precision 
of the extrapolation to zero frequency. On the other hand, the 
corrections due to the finite viscosity of the electron liquid do 
lead to a small but systematic decrease in the conductance, 
which is about 5% in the examined range of relative conduc- 
tances^'*. While this is smaller than the eiTor bar of the extrap- 
olated conductance, since the viscosity coiTection itself does 
not involve extrapolation (Eq. l24b . the resulting error bar of 
the relative correction, shown in the lower panel of Fig. |4l is 
noticeably smaller than the correction itself, and thus the cor- 
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FIG. 4: (Color-online) The Kohn-Sham conductance Gy with error 
bars resulting from uncertainty in the extrapolation to zero frequency 
(upper panel) and the viscosity-corrected conductance G^^ (shown 
relative to Gq^) (lower panel), as a function of the vacuum width d. 
Only the corrections due to the non-local kernel lead to changes in 
the conductance. The ALDA kernel does not affect the Kohn-Sham 
result within numerical error. 



dauer formula for broadening or opaque nanojunctions. Fur- 
thermore, we have utilized our formulation for examining the 
exchange and correlation effects on the conductance within 
the time-dependent density functional theory. We have shown 
that the long time limit determines the functional form of the 
exchange-correlation kernel that can lead to nonzero correc- 
tions. The exposed theory can be also used for ab-initio cal- 
culations for which the achievable precision in conductance 
calculation for a given size of a finite-size simulation cell has 
been demonstrated on a simple metal-vacuum-metal junction. 



Acknowledgments 

This research was supported by the Slovak grant agency 
VEGA (project No. 1/2020/05), the NATO Security Through 
Science Programme (EAP.RIG.981521) and the EU's 6th 
Framework Programme through the NANOQUANTA Net- 
work of Excellence (NMP4-CT-2004-500198). 



IX. APPENDIX: EVALUATION OF FUNCTIONALS ^ 
ANDff. 
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FIG. 5: (Color-online) The relative difference between the 4-point 
and 2-point Kohn-Sham conductances. The difference between the 
two diminishes for opaque barriers. 



rection is significant. 

Finally we numerically demonstrate the appropriateness of 
the 2-point conductance, G^^ = '^[x] for opaque barriers. The 
graph in Fig.|5]clearly shows that for d>4 the ^[ ]-dependent 
prefactor is very close to 1, which supports our arguments in 
Sec. [V] as well as is in agreement with previous numerical 
worb^ 



VIII. CONCLUSION 

In conclusion, we have presented a unified formalism, 
based on singular character of response functions, that gives 
the conductance of a general system of interacting electrons. 
We have derived a closed formula for the 4-point conductance 
in terms of the many-body response functions: the irreducible 
conductivity or the irreducible density response. The formula- 
tion allows for clear demonstration of the validity of the Lan- 



The functionals ^ [ ] and [ ] were given explicitly in the 
Fourier-transformed form in Eqs. (|9]l, ( fTOl i and Eq. O respec- 
tively. In numerical calculations it is more advantageous and 
numerically stable to use their representation in real space. To 
achieve this, we use the inverse transform of the irreducible 
density response function defined by 

Xiq,q')-J^e-'''xiz,z')e"''^'. (27) 
Several times we will have to resolve the integral of the type 

fizo) = [ %'^''>xiq,q')- (28) 

Since x{Qt) ^ 1 for small q, the integral is well defined and 
we can choose to interpret the apparent singularity l/q as 
l/{q + i5) or l/{q ~ i5) with 5^0+. Taking the former 
(the final result is independent of this choice) and using the 
inverse transform dZTb we find 

fizo) = ~ij dzdz'eizo~z)xiz,z')e'''''\ (29) 

where (z) is the unit step function. 

Using the definition of J# [ ] and using the integral ( |29] l twice 
with zo = we readily obtain 

^[X] = lim a dz [ dz!x{z,z!;ia). (30) 

a^0+ J-oo Jo 

The real-space form of the ^ is obtained using again the 
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Fourier transform: 



dq f dzdz' , M ijj 



qq' J 2n 



since we can Taylor expand e "> ' with the linear term giving 
the only nonzero contribution. 



q'=0 



on ^ 



= -ia / dz / dz!x{z,z!)—r- 
Jq J-oo q' 



a dz dz'x{z,z')z' 



(31) 



(32) 



(33) 
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For simplicity of presentation we assume that the induced field 
inside the electrode is absent, as is exactly true for jellium elec- 
trodes, used later in Sec. lVIII Generalization for an electrode with 
atomic structure is straightforward: The induced field E' is split 
into the local field induced inside the electrode, E'-'' (which has 
an average value zero over the electrode's unit cell), and the dif- 
ference = E' — E'-'' . The latter then defines the potential drop 
according to Eq. [3] and the former needs to be included into the 
response entering the functional ,^'^[d] where d = a -k { \ + x"" ) 
yNilhx' = 8E''<'/8E""\ 

From here on we use the atomic units with h = nip = e = I. 



